COSMOLOGICAL POST-NEWTONIAN EXPANSIONS TO ARBITRARY ORDER 



TODD A. OLIYNYK 

Abstract. We prove the existence of a large class of one parameter families of solutions to the Einstein- 
Euler equations that depend on the singular parameter e = ut/c (0 < e < eo)i where c is the speed of 
light, and vt is a typical speed of the gravitating fluid. These solutions are shown to exist on a common 
spacetime slab M = [0,T) X T^, and converge as e \ to a solution of the cosmological Poisson- 
Euler equations of Newtonian gravity. Moreover, we establish that these solutions can be expanded in 
the parameter e to any specified order with expansion coefficients that satisfy e-independent (nonlocal) 
symmetric hyperbolic equations. 



1. Introduction 

Einstein's general relativity is presently the most accurate theory of gravity. To completely determine 
the gravitational field, the Einstein field equations must be solved. These equations are extremely complex 
and outside of a small set of idealized situations, they are impossible to solve directly. However, to make 
physical predictions or understand physical phenomena, it is often enough to find approximate solutions 
that are governed by a simpler set of equations. The prime example of this is Newtonian gravity which 
approximates general relativity very well in regimes where the typical velocity of the gravitating matter 
is small compared to the speed of light. Indeed, Newtonian gravity successfully explains much of the 
behavior of our solar system and is a simpler theory of gravity that is less difficult to solve. By generalizing 
Newtonian gravity to the cosmological setting [25], it appears that the Newtonian theory can accurately 
describe gravity on all scales except in regions near compact neutron stars or black holes [9, 13]. 

Although Newtonian gravity is quite accurate, there are many situations where more accuracy is re- 
quired and general relativistic effects must be included. The post-Newtonian expansions were developed 
as a method to include relativistic corrections to Newtonian gravity such as gravitational lensing ef- 
fects, and energy loss through gravitational radiation. The goal of the post-Newtonian expansions is to 
approximate solutions to the Einstein field equations by a series expansion in the small parameter 

Vt 
c 

where c is the speed of light and vt is a typical speed associated with the gravitating matter. A number 
of formal calculational schemes for computing the post-Newtonian expansions have been developed and 
are widely used in both the asymptotically flat and cosmological settings [1,5,7,8,10,19,28,29]. These 
post-Newtonian computational schemes are one of the most important techniques in general relativity 
for calculating physical quantities for the purpose of comparing with experiment. Implicitly, they rely on 
the assumption that there actually exists solutions to the Einstein equations that admit post-Newtonian 
expansions to a certain order. This leaves open the question of existence of the post-Newtonian expan- 
sions. In view of the importance of the post-Newtonian expansions, this is a question of considerable 
interest. 

The main difficulty in determining when, and in what sense, the post-Newtonian expansions approx- 
imate a relativistic solution is that the region of validity for the post- Newtonian expansions is where 
e — Vt / c is close to zero. Therefore to understand the post-Newtonian expansions, solutions of general 
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relativity must be examined in the limit that e \ 0. In this limit the Einstein field equations become 
singular as the field equations contain terms of the form 1/e that become unbounded as e \ 0. The first 
person to overcome this difficulty and establish the existence of the 0**^ order post-Newtonian expansion 
(i.e. the Newtonian limit) without any special assumptions such as spherical symmetry was Rendall 
in [24]. In this paper, Rendall established the existence of a large class of one-parameter families of 
solutions to the Einstein- Vlasov equations that converge as e \ to solutions of the Poisson-Vlasov 
equations of Newtonian gravity. Using a different method based on techniques pioneered and devel- 
oped by Klainerman, Majda, Kreiss, and Schochet [14, 15, 26, 27] to study singular limits of first order 
symmetric hyperbolic systems, we established a similar result for the Einstein-Euler equations [20] and 
subsequently generalized it to prove the existence of solutions to the Einstein-Euler equations that admit 
post- Newtonian expansions to the 1*** order [21]. We also note that an improvement to at least the 2""^ 
post-Newtonian order is possible using the results of [22]. 

In this article, we adapt the methods of [20, 21] to the cosmological setting. We recall that the 
Einstein-Euler equations, which govern gravitating perfect fluids, are given by 

where 

T'^' = {p + c-'^pyv' + pg^^ , 
with p the fluid density, p the fluid pressure, the fluid four- velocity normalized by v^Vi = — c^, c the 
speed of light, and G the Newtonian gravitational constant. As shown in [20], these equations, upon 
suitable rescaling, can be written in the form 

(1.1) G'^ = 2e^T'^ and V,T*^' = 0, 
where 

(1.2) T^^ = {p + e^pWv' + pg'^ and v'v, = 

In this formulation, the fluid four- velocity v', fluid density p, fluid pressure p, the metric g^j, and the 
coordinates (a;*) i = 0, . . . , 3 are dimensionless. In this article, we restrict ourselves to cosmological 
spacetimes of the forrrj^ M — [0,T) x T^, and we will always use x = (x^) 1,2,3 to denote the 
standard periodic coordinates (with period 1) on the torus — x x . The coordinate a;° will 
parametrize the interval [0,r) and t = x'^/vt will denote an absolute Newtonian time coordinate. By a 
choice of units, we can and shall set = 1, in which case t = x^. 

The main aim of this article is to present a proof that establishes the existence of a large class of 
one-parameter family of solutions {gfj,vl, pe,Pe) < e < eo to the Einstein-Euler equations (jl.ip that 
(i) exist on a common spacetime slab M = x [0,T), and (ii) can be expanded in the parameter e to 
any flnite order £ e Z>o with expansion coefficients that are e- independent. For a precise version of this 
statement, see Theorem 11.11 below . To agree with standard terminology, we will say that these solutions 
admit a (i?/2)*'' order (cosmological) post- Newtonian expansion. 

In light of the significant and well-known difficulties that are encountered at both the formal and 
rigorous level in trying to develop post-Newtonian expansions on asymptotically flat spacetimes beyond 
the order 2.5 [23], it is somewhat surprising that these difficulties are absent in the cosmological setting. 
On asymptotically flat spacetimes, the problems that occur in the higher order post-Newtonian expansions 
are often attributed to the reaction of gravitational radiation with itself and matter. The analysis 
contained in this paper shows that this is not the complete story as the these effects are also present in 
the cosmological setting, but do not cause similar difficultiefl 

^The results of this article can be extended to other interesting spacetimes. A general discussion of spacetimes that 
admit Newtonian limits or post-Newtonian expansions will be presented in an upcoming paper. 
■^We will address this discrepancy more thoroughly in a separate article. 
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In the cosmological setting, the simplest one-parameter family of solutions that admits a post-Newtonian 
expansion to any ordei@ is the FLRW dust solution given by (^f.- = hij,v\ — = fJ.,Pe = 0) where 



(1.3) 
(1.4) 

(1.5) 
(1.6) 

and a(t) satisfies the equation 
(1.7) 



a'it) 



a"(t) 
a{t) 



ait) 



af{t)_ 
a{t) 



+ -ii(t) ^ 0. 



As is well known, the differential equation (II. 7p can be integrated explicitly to give 



(1.8) 



a{t) = ao ^ — 



i+\h^] ' and fx{t) = -[t + J — 



3^0 



where oq = a(0) and fiQ — fi{0) are positive constant^ 

In this article, we will, for simplicity, always assume an adiabatic equation of state of the form 



(1.9) 



P=fip), 



where / : R>o IR>o is an analytic function that satisfies f'{p) > for all p > 0. Even though the 
FLRW dust solution (|1.3p - (|1.8p does not arise from an equation of state of the form (|1.9[) . it plays the 
role of defining a Newtonian background that is essential for analyzing the limit e \ and generating the 
e-expansion. The role of the FLRW solution is easiest to see at the 0**^ order where it is used to define 
the cosmological Poisson-Euler equations 



(1.10) 
(1.11) 
(1.12) 



dtp + di[w p) 



3a'(t)o 
2 a{t) ^ 



p(dtW^ +w'diW^) = —(d-^f{p)+pd-^^ + a'{t)pw') 

a{t) 



A$ = a{t) p 



pd'^x 



(A = d'-'didj), 



of (cosmological) Newtonian gravity. We note that these equations agree with the Newton-Cartan field 
equations for a gravitating fluid formulated in adapted coordinates [17,25]. 

The cosmological Poisson Euler equations can be brought into a more familiar form by introducing 

^ 

Galilei coordinates [17, 25]. This is done as follows: suppose {p{t, x),w^ {t, x), $(i, x)} is a solution of the 
cosmological Poisson-Euler equations pTTUll - prT^ on M = [0,T) x T^. Then, letting M = [0,T) x 
denote the covering space, we define a diffeomorphism on M by 



: M — > M : {t,x) i — > {t,x/ ^/a{{}). 



■^This is an exceptional solution. Generically, it will only be possible to expand dynamical solutions to a fixed finite 
order in e. 

^With these choices, the big bang occurs at "time" t = 
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Lifting the cosmological Poisson-Euler equations to M, and then puUing back by ip shows thalQ 

(1.13) p{t,x)=p{t,x/^^), 

(1.14) wJ{t,x) = ^){h-'{t,x/^) + 



(1.15) <^{t,x) ^<i>{t,x/y^a{t)), 
satisfy 

(1.16) dtp = —w^ dip — pdiw^ , 

(1.17) dtw'' = -w^diw-' - -d''f{p) + g\ 

P 

(1.18) A4 = 4^p-y" pd^x 
where 



g-^ ^-^d-^i>~^\^l pd^x]x^ 



Defining a Newtonian potential by 

a short calculation shows that $ satisfies the Poisson equation 

(1.19) A$ = p 
while the acceleration due to gravity g"^ takes the familiar form 

(1.20) g'' = 

Together, equations ()1.16|) . p.l7|) . ()1.19|) . and ()1.20|) demonstrate that solutions to the cosmological 
Poisson-Euler equations determine solutions to the standard Poisson-Euler equations on the covering 
space M. 

To understand the limit e \ 0, we use a slight variation of the approach using in [20] and replace the 
metric gij and the fiuid velocity with new variables that are compatible with the limit e \ 0. The new 
gravitational variable u^^ is defined by 

(1.21) ""^TTO?! 

where 

(1.22) q'^ = h'-' + eV^ J/m", and J| = eSlS"", + 5\5l, 
while the new fluid four-velocity is defined by 

(1.23) = {I + ew'^)5l + 5\w^ . 

From these formulas it is not difficult to see that u*^ and are equivalent to the metric gij and the fluid 
four- velocity u*, respectively, for e > and are well defined at e = 0. 



'^In the Newton-Cartan theory, the fluid velocity 3-vector is the spatial part of a 4-vector w = dt + w^dj [25]. The 
formula 111.141 1 follows from the calculating the spatial components of ii = ip*w. The other two formulas 11.131 1 and 11.151 1 

- 

just follow from the definition of the puUback, i.e. p = ip* p and <& = i/j'^. 
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Following [20,21], these variables combined with a harmonic gauge can be used to cast the Einstein- 
Euler equations into a singular (nonlocal) symmetric hyperbolic system of the form 

A°{t, e, W)dtW = -C'djW + A'{t, e, W)djW + F{t, e, W). 
e 

This equation has the necessary structure that is required to use the results of [14, 15,26] to study the limit 
e \ 0, and to use Kreiss' bounded derivative principle [4, 16,27]to generate e-expansions. The beauty of 
Kreiss' bounded derivative principle is that it reduces the problem of generating order expansions in 
e to that of constructing initial data for which the time derivatives {dfW)\t=o are bounded as e \ for 
p = 1,2,...,^+1. This process of choosing such initial data is called initialization. Thus the existence 
problem for the post-Newtonian expansions is replaced by the problem of establishing the existence of 
initial data that can be properly initialized. The main technical result of this paper is to construct a large 
class of initial data for the Einstein-Euler equations that can be initialized to any order. The method 
for constructing this data is based on a technique introduced by Lottcrmoser in [18] who was the first 
person to prove the existence of a one-parameter family of initial data that depends analytically on e and 
converges to the expected Newtonian initial data as e \ 0. 

It is worthwhile to note that initializing the data to arbitrary order does not seem to be possible in the 
asymptotically flat case although it is possible for low orders [12,21]. In any case, initiahzing the data, 
whenever possible, provides a constructive method for generating initial data that has minimal initial 
gravitational radiation to an accuracy that is governed by the order of the e-expansion. This is certainly 
a significant improvement over other ad-hoc methods for choosing initial data that have been previously 
used in the literature. 



1.1. Notation. Before stating the main result of this article, we first introduce a number of function 
spaces. Given a finite dimensional vector space V, we let H''{V) denote the standard Sobolev space of 
valued maps on T'^. When = K, we just write H''. The only two vector spaces that will be used in 
this article are and the space of symmetric matrices Sat = { ) G Matxjv | u*"* = }. Letting 

(i/'i|'02)l2 = / ■il;itp2d^x= f f f tpi{x'^,x'^,x^)tp2{x'^,x^,x^)dx^dx^dx^ 
Jt^ Jo Jo Jo 

denote the standard inner product, we denote the projection operator onto the orthogonal com- 
plement of the constant function 1 by 

(1.24) n(V') = V- 1 y-ipeL^iT^). 

Given {ea}a=i any basis for V, we use the projection (ll.24p to define 

N 

H%V) = { V(a;^) = e H'{V) | (1|V'")l^ =0 for a = 1, 2, . . . , TV }. 

a=l 

We also define the standard hyperbolic evolution spaces 

Xt,ls{v) = n CP([o,r),i7^-p(T/)). 

p=0 

and write Xt,i,s if = K. Finally, given two Banach spaces X, Y and a open set U C X, we let C'^{U, Y) 
denote the set of analytic maps f : U ^ Y. 
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1.2. Main Theorem. The following Theorem contains the precise statement of the existence of post- 
Newtonian expansions to arbitrary orders. The proof of the Theorem can be found in Section [4.21 



Theorem 1.1. Suppose i G Z>o, s > 3/2 + i + 3, Cs is the constant appearing in the inequality (j2.77p . 
S - Mo/(2C,), y}^ e Si, y'/ £ S3, po e BsiH'), G H'{M.^), u'/ G H'+'-'iSs), u^, G H'-'iS^), and 
let To be the maximal time of existenc^ of solutions of the Poisson-Euler equations (|1.10p - p.l2p with 

initial data p\t=Q ^ Po + Po o,^d W"'^|f=o = WQ/(/io + Po)- Then for any T < Tq there exists an eg > 0, and 
maps 

p, G XT,e+2,s, K G ^T/+2,^(K'*) < e < eo, 

: k\ dik\ dM' e ^T/+2-p,s-p(§4) p - 1, 2, . . 
p G XT.e+2~p,s-p, e Xt,£+2-p,s-p(R'^) p = 1, 2, . . . £, 

: dtul^ G Xt,i,.-£-i(§4) G ^>i+i x (0,eo), 

e XT,i^s~e-i, w[ G XTa,5-£-i(K'') (g, e) G Z>f+i x (0, eo), 
AoGC"((-eo,eo),M), Zo' e C-((-eo, eo), Ao(0) = Zo'(O) = 0, 

such that 

(i) the triple Pe, w*} determines, via formulas (|1.2ip - (|1.23p . a unique solution to the Einstein- 
Euler equations p.ip m i/ie harmonic gauge for < e < eo on the spacetime region {t = ,x^) G 
M = [0,T) X T3, 

(ii) 

(l|u:^|t=o)L^ mul'\t=o)L^=yi\ 

n(ariuf^Uo) - .^u.^ii, n(ar^uf Uo) = ^^^^i, 

Pe|t=o = Mo + Ao(e) + /5o, aTid wf|t=o = (20(e) + Wo)/(Pek=o) 
/or < e < eo, 

(iii) {p, is the unique solution to the Poisson-Euler equations (I1.10p - (jl.l2p with initial data 

p\t=o = Po + P and w'\t=o = Wo/il^Q + P). 

(iv) forp = 1,2, ... ,1, {u^^ , p,w^} satisfies a linear (nonlocal) symmetric hyperbolic system that only 

depends on and {p, , $, u*-' , p, \ q ~ 1,2, . . . ,p — 1} , 

(v) forp G Z>j»_|_i, {u]P ,^p^,w\} satisfies a linear (nonlocal) symmetric hyperbolic system that only 

depends on e, {p, ,^,u^^ ,p,w^ \ Q = 1, 2, . . . , ^} and {u^^^ ,p^, w\ \ q = l^l,l^2,...,p — 1^, 



^See Proposition 14. II 
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(vi) {ul^ , p^,wl} and {ul^ ,p^,wl} (p > i + 1) satisfy the estimates: 

\\unt)\\Hs+.+e\\d,unt)\\H^+e\\mul'{t)\\H.-.+e'^^^^^ < 1, 

Wp^mns + HmH^ + WdtpMlH^-^ + \\dtwl{t)\\Hs-i < 1, 

fpMiH^-^-^ + \\wi{t)\\H.-e-^ + wdtkmH^-'-^ + iia/«:(t)iiH=— < 1, 

for all {t,e) G [0,To) x (0, eo), and 

(vii) {iT'^ , Pcwl} admits convergent expansions (uniform for < e < eo) of the form 



^ °° 



p=l 







p=i+i 



°° 



p=i 



d'(p, = dfp + Y,ePdfp+ E ^"^t'^e, 



P=i p=^+i 
where the expansions are convergent in C*^([0, Tq); i?"^*"^^^"^) for = 0, 1. 
An important point not explicitly stated in the above Theorem but follows from the proof is that for 



p= 1,2,. 



the equations satisfied by coefficients {u*-', p, w'} from Theorem ll.il fiv) can be derived by 



assuming a harmonic gauge (j2.18p and substituting the expansions of Theorem 11.11 (vii) in the Einstcin- 
Euler equations (|l.ip and collecting terms together of the same order in e up to order i. It is this fact 
that guarantees that the expansions of Theorem 11.11 (vii) coincide with the post-Newtonian expansions 
of order 1/2. 

As a final remark, all of the results of this article can be adapted to include a cosmological constant 
A. The basic changes needed to be made include replacing the stress energy tensor (11.21) by 

- (p + e^p)v'v^ + pq'' - ^g'\ 
and replacing the FLRW equations (|1.5p and (|1.7p by 

'm ~ 

a{t) 



p{t)+A=- 



and 



a{t) 



a'(t) 



ait) 



+ -(M(t)+A) =2A. 



As shown in [3], the inclusion of a positive cosmological constant guarantees the long time existence 

Q Q 

of small perturbations of the constant density solution {p ~ p,w^ — 0,^ — 0} of the cosmological 
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Poisson-Euler equations (I1.10p - (ll.l2p . The importance of the long time existence of solutions is that it 
is a necessary ingredient of any analysis of a lower bound on the time of validity of the post-Newtonian 
expansions as a function of e. We plan to address this problem of determining a lower bound on the time 
of existence of the post-Newtonian expansions in the near future. 

2. The Einstein-Euler equations 

2.1. Reduced Einstein Equations. In order to derive a suitable symmetric hyperbolic system for the 
gravitational field equations, we introduce new coordinates related to old ones by the rescaling 

and let 

In the new coordinates, the spacetime metric g.ij and the FLRW metric hij (see (jl.3p ) are given by 



-g,^ = JW^9^, and = + a{t)5i5i5 



IJ, 



respectively, where Jf is defined in (11.221) . For latter use, we record the non-zero independent components 
of the Christoffel symbols 7^^ and the curvature iZijki of the metric hij: 



(2-1) - 2 ait) ' 

(2.2) 7?/ - 



(2.3) 7^0/0/ — 

and 



2a{t)a"{t) - [a'jt)]^ 



(2.4) :^2121 = 7^1313 = 7^2323 = j[a'(0]'- 

As discussed in the introduction, we take the symmetric 2-tensor u^^ as our primary gravitational 
variable where 

(2.5) = ^ (^^^-9'' - ^''^ ^ VWla"' - V^(^*' + e'w*'), 

and 

(2.6) I5I -det(.g,y) and \h\ ^ - det(h,j) = [a{t)f . 
Observe that the metric can be recovered from the u*^ by the formula 

v/~|/i|det(g'=0 

where 

(2.8) g'^ = h'^ + e'^u'^ . 
Substituting (|2.7p in to the standard formula for the Christoffel symbols gives 

(2.9) f ^- = 7^^ + e2 + iff'='5^m5jn5zS™" - iff^'s^jSmn^i^"" + ^gi„,Sf^D ^)u'"') , 
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where (cjij) ~ (g^^) ^ and Dk is the hij covariant derivative. Using this formula, the Einstein tensor G*-' 
of the metric gij is given by 



(2.10) 


— 1 /^4 7 


4n 


g^^DkDiu'i + e2(a^ + 4^ + 4) + h'^ + ed^ + e'^i + Ae^f'' 


where 








(2.11) 


ij 

a{ = 


\{\9kl9-mn - gkmgin) {o'" ~ \g'= 9^-^) D pu'^' D ,U^^ , 


(2.12) 


ij 
0-2 = 




(2.13) 


ij 
o 




(2.14) 




g'WkDiu'^' -2DiDku''^'g^^\ 


(2.15) 


— 




(2.16) 


ij - 
^2 — 


1 9 

— ij —kl n~) \ Tf (i —j)k —hn 


and 








(2.17) 


fJ = 






To fix the 


gauge, 


we set 




(2.18) 









For e > 0, it is clear from (|2.5p that this is equivalent to 

and this is easily seen to be equivalent to the harmonic coordinate condition 

Defining the reduced Einstein tensor G^^ by 

(2.19) G% = - = hig^'OkDiu'i + e^{a\' + 4 + 4) + e^| + e'^i + Ae^r^), 

e \n\ 

the Einstein equation G'^ = 2e^T'*^ in the gauge (|2.18p becomes 

(2.20) G)j=2e2l£|r^ 
where 

(2.21) r^' = (p + e2p)i;«^^^' and w^w^ " 

To write the reduced Einstein equations p.20p in first order form, we introduce the variables 

(2.22) u'^ =eu*J' and u!^ ^ DkU^K 
With these variables, we have that 

or equivalently 

(2.23) d,u!^ = + 27L4)'" - 27L^.f - -nuJ^u^^--. 
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In particular, this implies that 

(2.24) Bon"^ = Bjv^ + 27L"o^'" " 27^"^" " -noiJ'u^^"-, 
and hence 

^k^iu — 9 OoUq +9 OqUj +9 OjUq +g OiUj+g (-^ ^lU^^ + ^^u^ ) 

(2.25) = f°doul' + 2eu°'diui,' + g" Biu'j + ed^{ + e'cti , 
where 

(2.26) d-^ ^h'^(-lj^u^^ + lj<^^up^ 



e e 



and 

(2.27) 4 = (-^7S< + 'CW' ") + ("tL^^o^" - 'iL-r - 1^0/™^'^^")") . 
Setting 

(2.28) ey = -7j™< - -7L"j " ^^ojm^ , 

equations (|2.24p and (|2.25p can be used to write the reduced Einstein equations (|2.20p in the following 
first order form 

-a{t)g°°Btul' = 2a{t)u"'aiul' + ^a{t)g"div'] + a{t){^^ + d\') 

(2.29) + eait) {a\' + 4 + 4 + ^i + 4) + 4ea(t) [t^^ - Mr^'^ , 

(2.30) am^dtv'j = -^a{t)g'''djul= + a{t)g"e% 
and 

(2.31) dtu^' ^ul' --^^gu^^''. 

For their definition, the reduced Einstein equations (|2.29p - (|2.3ip require that the matrix g*-' is invert- 
ible. For fixed — •\/2/(3/xo) < tq < and ri > 0, it is clear from p.Sp that 

< a{To) < a{t) < aiTo) Vte[ro,Ti]. 

This implies that the set 

Vro,n = { {r'n e M4X4 I det{h'^ + r'^) > V < G [tq, n] } 

is open and contains the origin (r*-' ) = 0, and moreover, that the reduced Einstein equations (|2.29p - (|2.3ip 
are well defined for all t S (tq, n) and (eu*-') e Vto,ti- 

2.2. Regularized Euler equations. In the coordinates (a;*), the Euler equations are given by 

(2.32) Vjf = 

where T*^ = (p + e^p)v^v^ + pg^^ and the fluid velocity is normalized according to 

(2.33) v,v' = ~\. 

To derive a symmetric hyperbolic system for the Euler system, wc follow the method of [2] and differentiate 
((233)) to get 

(2.34) w.Vjw' = 
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which in turn implies 

(2.35) v^ViVjv' = 0. 
Writing out (|2.32p explicitly, we have 

(2.36) (B^p + e^d^p)v'v' + {p + e^p){v^V^v' + v'V.v^) + g'^d^p = . 
Next, we observe that the operator 

projects into subspace orthogonal to the fluid velocity i.e. l-i\L\ — L-j, and Ljv^ — 0. Applying this 
operator to project (j2.36p into components parallel and orthogonal to v' yields, after using the relations 
([233)1 - ([235)1 ■ the following systenfl 

(2.37) .^M^y%p + f'(Pj L)V,v^ = , 

{p + e^f{p)y p + e^f[p) 

(2.38) M^.v'^Vkv' + -^^L'Ap = , 
where 

As discussed in the introduction, we need to introduce a new fluid four-vector by 

(2.39) = v' -e =v' - 

e 

So, letting 

(2.40) Mv^{p,w'f 
allows us to write the system p.37p - p.38p as 

(2.41) Aljdtv^^Ai,diw + FM, 
where 

/ f'{p){l+ew") ef'jp) jQ \ 

(2.42) Alt= (P+f^nP)y P+^^fiP^^ , 

/ f'ipW , f'ip) 7-j\ 

(2.43) L'Tirrr ""3^^ ' 



and 



(2 44) Fm = i P+'^^^P>^^^''^ - P+^''f{p)^J^a^ 



^Recall that we are assuming that the fluid satisfies an adiabatic equation of state p = /(p) (see II1.9II ). 
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Next, a straightforward calculation using (12. 7p and (|2.9p shows that 
My = % + 2,50(50 + em,, (e, t, u'^^w'^), 



Hi 



2 a(t) 



1 



a{t) 



where the maps rn^j, t^, q^, qi, and q^ are analytic in all their variables provided that t G (to,ti) and 

(eu*^) G Vt,ti- Using these expressions, we can decompose A\,j, A{,j, and F]\f as 

(2.45) Al, = AO, + eA°i^,{t, e, i.'^ p, u;^, 

(2.46) Ai, - + eAij_,{t, e, p, u;^), 

(2.47) Fm = Fm,o + eFM,i(i, e, Ufc*', P, w''), 

where 
(2.48) 



(2.49) 

(2.50) 

Fm,o 



fip) 




2^o^o 



f'(p)3a'{t) 
p 2 a{t) 

{h, + 2S'^S]) + ij^0(3C _ a(t)5KL^.^^) - kS^' (^uf + S^^uf^) - 



and the maps Ajjj q, i^M,i are analytic in all their variables provided that t G (tq, ri) and (eu*-') G 

2.3. A nonlocal symmetric hyperbolic formulation. To bring the reduced Einstein equations (|2.29p - 
(|2.3ip into a form that is suitable to analyze the limit e \ 0, we replace the uj with the variables 



(2.51) 

where the ^^-^ satisfy 

(2.52) {6'' + a{t)eu'-')didj^'^ = 4e'a{t)n I - 

\vA 

A short calculation shows that 



\9\, 



(2.53) 4e2 ( _r«i ]=^{^p~p,(f))8l8l + ^{%b'ipw^^ ^\8l8lphuu^'\ + e^S}^ {e,t,u^\ p,w% 
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where again the map Sq is analytic in all variables provided that t e (to,To) and (eu*^) G Vto,ti- 
addition to , we will also need the time derivative 



In 



which satisfies 
(2.54) 



\h\ 



\9\. 



\h\ 



Using (|2.37p - (|2.38p to replace the time derivatives of p and w"^ in favor of spatial derivatives, we find that 



\9\, 



\h\- 



3 a'{t) 



^^ dt - T'M = - 4 di{pw^) + -^{p - Pit)) SX + e 



2 a{t) 



riw'S^^pw^^] 



ait) 



Si'S^^^djfip) 8^ (^4V«^^') + si'S^^pw'^ + %p5U'^ 
'diipw') - \^P) hiu'^ + P {a'it)Sjju''' + hi4') 



(2.55) 
where 



SI' (e, t, u''\p, w'') + SI' (e, t, u'^\p, w^OkP, Okw^u'^!,) 



P' = -u^ + \S''{H" - ait)SKLU^'^) ~ IS'' (^^^.?° + SKLuf^^ , 



and the maps S*^^ (a = 1,2) are analytic in all variables provided that t e (To,ro), (eu*^) G Vtq.Tq, and 
1 + > 0, and S^' is linear in (S^p, Okw'^ , u'^\) . 
Substituting ^jl^ into ([2:291) - ([230]) gives 



(2.56) 
(2.57) 
where 
(2.58) 



iit)g°°dtul' = 2ait)u"diul' + ^ait)g''' diWy + ait){c^ + d^) 

+ eait) ia\' + + ^ + 4 + 4) + , 



ait)g''dtWy = -ait)g''djul' +ait)g 



(j)''it) = 4ea(t) ^1 



rj-ij _ ffl rpij 
\h\ 



L2 



Differentiating (|2.58p with repsect to t while using ()2.53|) and (12. 55^ . we find that 



(2.59) 
where 



<l>'"{t) 



-6^^it) + Fi' ~ell 



a!if)-Si^ +ait){S^ +S^) 



L2 



--^ait){l\pw')L^5'^i5\^ - 8a(i)^l p (^-S^^u'^'' + i4'5^')"(3C - a(i)<5KLU^^) 

(2.60) + 2(np|4'<5J")^(u?" + ait)6KLuf''))L^ - ^a{t){l\p{a' it)6iju" + hkml')) l-SIS^ 
Next, we define 

(2.61) w = (M;^ly;^u'^0'^w)^. 

Then it follows from ((23T1) . (jTiTj) . ((2:45)) - ([CTf)) . ([236l) - ((237|) . and ({2391) that W satisfies 

(2.62) A°atVl^ = -C^diW + (A^ + ev4()a7W^ + i^o + e-F'i, 
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4° - 

A.Q — 



^M,0 + ^^M,l 

/a(t)(l-euOO) 







0\ 



(5" + ea{t)u") 
1 

1/ 



^0 



Ah 



(2a{t)' 



,01 



a{t)u" 







Am A 


































0/ 



ch 









5" 0\ 





0/ 



M,0 



Fm,o 

P 2 a(t)_ 



4ait)A-'IL{p-^^it)), 



{p - Pit)) 



Fi = 



a{t){aY + 4 + + c^' + d^i) 
7 (-Fm,o — Fm,o) + Fm,i 



- 1 



1,2 



^1] 
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2.4. Well-posedness of the nonlocal system. With the evolution equations in a suitable form, we 
now verify that the system is well-posed. To do this, we will repeatedly use the following elementary 
facts concerning analytic maps: 

Lemma 2.1. Let X , Y , and Z he Banach spaces with U C X and V C Y open. 

(i) If L : X — > Y is a continuous linear map, then L g (X,Y). 

(ii) If B : X X Y — > Z is a continuous bilinear map, then B S C"(X x Y, Z). 

(iii) Iffe C"(C/,r), g G C"^iV,Z) andr&nif) C V, then g o f e C^{U,Z). 

We also recall the well-known Multiplication Lemma. 

Lemma 2.2. Suppose si, S2 > > Q and S3 < si -I- S2 — 3/2. Then there exists a constant C > such 
that 

UMh^^ <CUi\\h^A\4^2\\h'2 

for all V-'i e H""^ and ip2 G H''^ ■ 

This lemma shows that the bilinear map 

(2.76) H'' X H''' 9 (V-i, V'2) ^ G H"^ 

is continuous, and hence analytic, provided that si, S2 > S3 > and S3 < Si + S2 — 3/2. In particular, 
iJ* is a Banach algebra for s > 3/2 which implies that there exists a constant Cs > such that 

(2.77) ||^/'i^2||. <C,||^i||h.||7/'2||h= 

for all V'l I V'2 G ■ This can be used to prove the following important proposition concerning analytic 
maps. For a proof, see Proposition 3.6 of [11]. 

Proposition 2.3. Suppose s > 3/2, F G C^(_B7^(]R^), M), Cg is the constant from the inequality (|2.77p . 
and that 

F{y,, . . . , yAT) - Fo + ^ c„ y'^^y^ ■■■y%" 

\a\>l 

is the powerseries expansion for F{y) about 0. Then the map 

{BR{H')f 3 (V^i, V2, . . . , V'at) ^ F{iP,,iP2, ...An)^H' 
IS m C'^{{Bn/cSH'')f,H'^ 

N 



\a\>l 



for all (^i,...,^jv) e (Bn/cAHn)''- 

The first step in establishing well-posedness is to show that the maps $ and $ are well-defined and 
analytic. 

Lemma 2.4. Suppose R> and s > 3/2. Then there exists an eo > and an analytic map 

(-60, 60) X (to. To) xH'x H^R") x Bn{H%S^)) 3 {e,t, p,w\u'^) ^ ($^^) G H'+^S^) 
that satisfies (|2.52p and 

^''l^o-^a{t)5l,siA-'n{p-^,{t)). 
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Proof. First we observe that for a fixed i? > 0, the Born series 



1 + ea{t)A--^u"didj] A^^ = e"(-l)"a(i)"(A-^u"9/9,/)"A- 



the Multiphcation Lemma [2.21 and the invertibihty of the Laplacian A : show that there 

exists an Co > such that the map 

(2.78) 

(-eo,eo)x(To,To)xBfl,(ff'+'(S3))xF^(§4) 9 (e, t, u", *^^) ^ [{S''' +ea{t)u")didjY' "f^' & H'+^S^) 

is well defined and analytic. Also, by Lemma 12.11 and Proposition 12. 3[ we sec that (shrinking eq if 

neccessary) the map 

(2.79) 

(-60, eo) X (ro,To) x x H%R^) x i7^(§4) 3 {e,t, p,w\u'') ^ 4e^a{t)U (Mf^^ - e iJ-(§4) 

is well-defined and analytic. The proof then follows directly from composing the two analytic maps 

mrn-nrm. i.e. 

-1 



h\ 



which is again analytic by Lemma [271] 
Lemma 2.5. Suppose R> and s > 3/2. Then there exists an eo > and an analytic map 
(-60, eo) X (ro,To) x x H'iR") x B^H^Si)) x H^S^) 3 {e,t, p,w\u'\u]^) ^ G H'+^S^) 

that satisfies (I2.54p and 



□ 



l a'jt) 
2 a{t) 

Proof. The proof follows from a routine adaptation of the proof of Lemma 12.41 □ 

Next, we introduce the space 

Hr = ff'(§4) X {H'{§i)) X BR{H'{§,i)) X §4 X iJ" X Br{H') x H'iM.^) 
and let W = ni^. 

Lemma 2.6. Suppose R> and s > 3/2. Then there exists an eo > such that the maps 

(to. To) X (-60, eo) xHj,3 {t,ey^ ,Wy p,w° ,w'f ^ F„ e {a = 0,1) 

are analytic. 

Proof. Follows directly from Proposition l2.3l and Lemmas 12.11 12. 4 [ 12.51 □ 

With the analyticity of the maps Fq,Fi established, local existence, uniqueness, and continuation of 
solutions to the nonlocal symmetric hyperbolic system (|2.61l) follow from standard arguments (see for 
example [30], Chapter 16). In particular, we can apply the local existence results of Schochet [26,27] (see 
also [14, 15]) to obtain the existence of solution to (|2.6ip on spacetime regions of the form D = [0, T) x 
where T is independent of 6. This will be discussed in detail in Section [ 
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3. Initialization 

3.1. Constraint equations. In order to solve the initial value problem for the Einstein equations, we 
must first construct initial data that satisfies the following constraint equations on the initial hypersurface 
defined by t = 0: 

(3.1) (G"^' - 2e'^f"^)\^^^ = (gravitational constraints), 

(3.2) L'iu'-' I = (harmonic gauge condition). 



and 



(3.3) 



(fluid 4- velocity normalization). 



A short calculation using (|2.ip - (|2.2p shows that the harmonic condition (|3.2p is equivalent to 



(3.4) 
(3.5) 



dtu 



00 



2a{t) 



{Su"^" + a{t)Suu"), 



e 2a{t) 



Using formulas (|2.7p and (|2.10p - (|2.17p . it is not difficult to verify that the gravitational constraint 
equations (j3.ip do not involve second order time derivatives. In fact, using p.4p - (|3.5p . the top derivative 
terms can be expanded as 



1 



a{t) 



'ju 



IJ 



(3.6) 



a' it): 



ait) 



a'it) 



-10 



2ait) 



where /o and /i are analytic in all of their variables, and /2 is linear in idxdLU^^ , Oku'^-' , dtu^^)- Together, 
(|2.7p . (|2.10p - (|2.17p . (|2.53p . and p.6p show that the constraint equations (when evaluated at i = 0) can 
be written as 



4(p - /i) + didju 



IJ 



o"o 



eao 



-Apiw^'Si + w^)+ 5'j ( didtu'-^ + 



00 



ao 



-/o 



2ao 



fiie,u'\dKdLu'^,dKu'\dtu'') 



(3.7) 



+ fi (e, v^^dtV^^dKV^^) + fi (e, u'\p, w\pw') 



= 0, 



where for any i? > there exists an ep > such that the maps fa (a = 2,3,4) are analytic in all their 
variables provided |e| < eg, and \u^^\ < R, /| is linear in idxdhU^^ ^Qku^-' ,dtu^^), and /| is quadratic in 
(9tu*-' , m'^ ). Also, an easy calculation using (|2.7p shows that (|3.3p takes the form 



(3.8) 



0,1 5ojW^ + y/e^igojw'^)^ - gooie^gijw^w'^ + 1) 



e.9oo 



where the map /o is analytic provided |e| < eo, | < R, and \w^\ < R. 
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3.2. / J-components of the reduced Einstein Equations. The /J-components of the FLRW wave 
operator acting on u^^ are given by 

1 



(3.9) 



2 a(t) 



'a{tf 



dt 



a!(t)\ a!^ 
a(t) J a{ty 



-J J 



2a(t)3 



A calculation involving the harmonic conditions ([Ojl - dSTSl) . and formulas (l2Jl) . (|2TT|) - (|2T7|) . (12719]) . 
(|2.53p . p.9p shows that the / J-components of the reduced Einstein equations can be written as 



Au" + edi 



a{t) 

,2 



ait) 



,0/ 



7 a'{t) 
2 a{t) 



dtu" 



a{t)' 



(3.10) 



where for any R > there exist an > such that the maps Pa {a = 0,1,2,3) are analytic 
in all variables provided —To < i < ri, |e| < eg, and \u^-'\ < R. Furthermore, p^"^ is linear in 
(e'^d'^u^^ , edtdMU^^ , QmQnu^^^ , dMU^\ (^dtu^^ , eu^^), p{-^ is quadratic in {dRU^-' ,(-dtU^^ ,£u^^), and P2'^ 
is linear in (m*-' ), 

3.3. The Euler equations. Directly from equation (|2.4ip and formulas (|2.45p - (|2.50p . it follows that 



(3.11) 



(3.12) 



dtw' 



dtp 



a{t)p ^ a{t) 4 



a{t) 



KL 



qi {t, e, p, , u''')+q\ (t, e, p, , u^'' , dip, diw^ , dtu^*' , dju^^) 



J 3 a'{t) 



90 



where for any i? > there exists an eq > such that the maps qa, q^ = 0, 1) are analytic in all variables 
provided tq < t < ti, |e| < eq, I^'-'I < R, and \w"\ < R, and qi,q{ are hnear in {dip,diw^ ,dtu^'^,diu^''). 

3.4. Higher order time derivatives. As discussed in the introduction, Kreiss's bounded derivative 
principle requires us to calculate higher order time derivatives of p, and w\ The fact that the 
constraint equations must be satisfied complicates this task, and we find it advantageous to introduce 
the following rescaled variables: 

(3.13) dfi^^u"^yj/ + e'ul-\ 9f = + 9,^=o^°° = « + > 0) 

p|t=o = A*o + Ao + /5o, w°\t=o^WQ, {pw^)\t=o = Zq 



(3.14) 
(3.15) 

where y^"^ , Zq, Aq are constants, and 



d: 



and d\ 



(3.16) 

We also define 

(3.17) 

(3.18) 



T3 



Pod X 



T3 



Wq (fix 



X3 



~0j j3 
Wq a X ■ 



w 



Wn 



> 1) 



ul-'cPx = (^>0). 



T3 



u/ — (Uq"' , m/ , . . . , u/), — (2/2 i---;2/£+2)i 



P\+1 = (Pl,P2, 



and 



^^+1 = (w1>2:--->^+i)- 
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Differentiating the harmonic conditions ()3.4p - ()3.5|) . the IJ components of the reduced Einstein equa- 
tions ()3.10|) . and the Euler equations (|3.11|) - (|3.12|) with respect to time and evaluating at t = while 
using the variables (|3.13p - (|3.15p yields the following system of equations 

(3.19) A4° - 4ao(po + Ao) - eFO(e, .yi' , u^, d^+'u'^ , Xo, Po, 4, w'o) = «, 

(3.20) AuSJ^ - 4ao(zo^ + w'^) - eP' {e, y^ ,y{' , ui^ Aq, po, zi, 4) = 0, 
^ui' - aoyiU - PoM'^ y^+\,dKnT:i^T: Ao, Po, h, 4: wf ) 

(3.21) - e'Pl^,{e,yl' ,y{\yf_^%,d^u\' ,u%,,d,uf+\,dMdNuf'^,uf^^^^ = 0, 



(3.22) 



,-,00 



E 

p=0 

1 dP+^a 



3dP+Mn(a) 
2 dtP+i 



SpoVo 



00\ 



Sij{y'e^p + e^ul:Lp) 



= 0, 



(3.23) ufl^+d.ul'^ + eY^f'Y- 

p=0 



5 dP+^ ln(a) 



t=o 



+ Spoy'o') = 0, 



(3.24) wo - eQ°{e, y'^' , ui\z^, w^) = 0, 

wl+i - Qo,^(Ao,a^po,9"p^,Zo,9"w;^a^w^) 

(3.25) - eQi^,(e, ,y{\ y',i,y/^,, Oku^ , Xo, d^po, d^Pt, 4, d'^w^dKuf) - 0, 

and 

p\+i - Qoj{Xo, d^po,d^Pi, z^, d^w^,d^w^^) 

(3.26) - eQl^£(e,yo^y",y"l,u^+l,9i^u^'', Ao,5^Po,9"p^,Zo,9"w^) = 0, 

where < |q;| < 1. Next, we set 

(3.27) vl^oo = ^^00 _ 4„^~^ _ ^nFO(e, y^,yi', utd^+'u^^Xo, po, 4, 4) 

(3.28) M'"'^ = AuS]^ - 4aozDo^ - ellF-' {e, y]^ ,y[' , d^+'u^^ , Xo, po, 4, ^o)' 
^i' = Au'/ - UP^j (j/OO, yf/i , 9;,ur , uT, Ao, Po, p„ Zq^, ri;o^, wf ) - 

(3.29) einP/;|(e, 2/^^ yf ^, y,^i ^i^w.^^, ^M^jvuf ^, uf+i Ao, po, P„ ^o"", , ^f), 



(3.30) TO" = Ao + -^{l\F°{e, yl' ,y^\u[^dKdLu'' ,dKu'^ , Xo, po, z^ w^)) 

4ao 

(3.31) Tl-' = zi + ^{\\P\e, 2/f ^, ui^ ^^aiit^^ 9^"'^ Ao, po, ^o^ .^o'))' 

4ao 

JL ^+2 - 2 ^ \1Ko,£(,2/o ,y£-l-l,CifUc ,Ao,Po,P<;,^o >^0 j/- 

ao 

(3.32) -(l|F/i(e, ^ yf+i, uf , a^u^^ tJ^iif/i, ^M^^i^f ^ <i Ao, PO, p„ ^o"", <>f )), 
ao 
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^^00 _ -00 



p=0 










2 dtP+^ 



t=0 



(3.33) 
(3.34) 



OJ -OJ 



^£+1 — 



p=0 



s\ 5 dP+i ln(a) 



2 dtP+^ 



'iiV-p + ^poyo 



t=o 



0J\ 



(3.35) 1]° = zI;o-eQSJ(e,yo'^^,uf ,2o'^,?I'o^), 

(3.36) - eQl,,(e, y[i^,n%„dKn\\\^, d^po, d^Pe, 4, d^w^ .dKnf) , 

^£+1 = p^+i " Qo,f (Ao, a^po, zq, a^w^ a^w^) 

(3.37) - eQi,f (e, y]^ ,v[\ y^:,^!, u^-'+i, ^xuf , Aq, 9"po, i9"p^, Zq: ^^w^O- 
Gathering all of the maps p.27p - p.37p together, we define 

(3.38) = 2, 1]^, n^,^,, TO", T0•^ T,^^^)'', 
with 

(3.39) *f = (vI/;^...,vI,;^), T[i, = {T'^\...Xi2l 

(3.40) = (f7j,...,l]^,^J, and = (fli, . . . , f),+i). 
We also define 

^ksA = Bni^i) X §3 X i?^+i-^(§3) X X B5(i?^) X Bfl(iJ^(R3)), 

£+2 £+2 
Br{W+'^) X X Sfl(iJ^+2(R3)) X i?^+2-P(R3) X Bb{H'+^ (§3)) , 



p=l 



£+2 



= 1 

£+2 



X Y[ H'+^-PiSs) X BRiH") X Jl H'-Pim*) X J| i/^^P x {-5,6) x x (§3^ 

p—i p—i p—i 



and 



£+2 £+2 

Zl = i?" X [I i7«+2-P X F"(R3) X J]^ i?"+2-P(M3) 

p— 1 p—\ 



-1 



^+1 



X Jl H'-P{S3) X i/^'+i X [| H'~P{R'^) X [| X R X X (§3)^ 

p—O p—1 p—1 



Proposition 3.1. Suppose i e Z>o, s > 3/2 + ^, i? > 0, is the constant from the inequality (I2.77p . 
(5 = fio/{2Cs) and set 

Qi ^ {y'o ,y" ,ulli,iill2, Po,Wo)'^ , and = {iii^, ^4+2, ^u)^, i4+2' Pi+2: ^0, ^^Y^U) ■ 
Then there exists an eo > such that the map 

i~eo,eo) X Xli s,£ x 3^fl,5/ 3 (e,6'<?,?7^) ^ e Z| 

is analytic. 
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Proof. By definition, 



p ^o + Ao + po Mo (1 + (Ao +/5o)/Aio) ' 
Since the map k{r) — 1/(1 + r/fio) is in C"^((— /io, Mo), K), it follows from Proposition 12.31 that the map 



- 

P 



is well defined and analytic for 8 ^ [i^j (2Cs). With this map well defined, we can recover from Zq and 
w\ by the defining relation = (zq + w}^l p. The rest of the proof now follows from a straightforward 
application of Proposition 12.31 and Lemmas 12.11 and 12.21 □ 

3.5. Initialization to an arbitrary order. We are now ready to prove that there exists a large class 
of initial data that can be initialized to an arbitrary order. 

Theorem 3.2. Suppose £ ^ Z>o, s > 3/2 + £, R > 0, and S = po/{2Cs) ■ Then for any 6^ ^ g ^, there 
exists an open neighborhood (— eg, eg) x Ug^ C M x s t '^f (0' (^o > 0, R > R) and a map 



m 



r.oo 



u, 



0.7 



'£+2: "^+2 7 



that satisfies 

Ee{e,ef,,r,e{e,ef,)) =0 and iB°(0, 0,) = A(0, 0,) = Zo'(0, = 
for all {e,0e) e (-eo,eo) x Ug^. 

Proof. We first establish that given a di £ Xf, ^ ^, the equation |c=o = has a solution. 
Lemma 3.3. For any 9e G X^ g ^, there exists a R > and a fje G 3^^/2 o e ^^^^ satisfies 

Proof. To begin, we consider the fixed data 



Next, we note that ^| 



e=0 



reduces to (see 



oo,s,e- 



(3.41) 
(3.42) 
(3.43) 
(3.44) 
(3.45) 

(3.46) 
(3.47) 

(3.48) 
(3.49) 
(3.50) 
(3.51) 



Ao = 0, 
0, 



-0^ 



4aow)o , 



Aw! 



-OJ 



Au^-' = IlP',-;^{{p,,w^, p„ 



yp+2 



f.OO 



"p+1 



ao 



9=0 





■3d«+Mn(a) 





2 



t=0 



00\ 



1 



2 



t=0 



Wq = 0, 

'Wp+i = Qap(Ao,Mo,Pp,2 



^ ,7,0 ,7,^ wJ r,oo 



"p 7' 



Mp+1 



= Qa^p{X(),Po,Pp,z^,w°,WQ,w^p). 
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where the maps 

Po,p : Bs{H') X H^iR^) x K x (§3)^"^ x H'+^ 
p p p 

X Y[ H'+^-'i X Yl H'-" X Yl H'-^R^) — > H'-P{S3) {0 < p < £), 

q=l q=l q=l 

P 

Qip : (-(5,(5) X BsiH') x Y[H' " xR^ xH' x H'{R^) 

q=l 

P P 

X Y[ H'-^iR*) X X Yl H'+'^-^iR*) — > H'-P{R^) (0 < p < ^ + 1), 

q=l q=l 
P 

Qo,p : (-<5,'5) X BsiH'') x J| H"-'' x R^ x x H'^iR^) 

q=l 

p 

xY[H'-'^{R'^) — rH'-P-^ (0<p<^+l) 

q=l 

are all analytic for 5 — /io/(2Cs). 

The invertibility of the Laplacian A : then implies, by equations p.43p - p.44p . that 

(3.52) = 4aoA-ipo and u^-' = 4aoA"^u'i(. 
Substituting these into (|05l) . dOe]) . ([QT]) . ([S^O]) . and ((33T|) (for (p = 0)) gives 

(3.53) u',' = A-inP5;'o(po,<,2/o^"o°), 

(3.54) y,^''^l(l|P^;^o(po,^o^,yo°^"^°)), 

ao 

(3.55) = -diul^ - ^(3u™ + 32/°'' + a(0)(5,,yo"), 

(3.56) < = Q^^o(Ao,po,4>o:"*'°)> 
and 

(3.57) pi = (3o,o(Ao,/5o,Zo,Wo). 
From ((OT)) and l|333)) . we then obtain 

(3.58) u?^ = -^/it^'^ 
Substituting into dSHD, we find 

(3.59) = A'in(po, ^'o'', 2/2^^ Pi>f ) • 

With the base case covered, we proceed by induction. So assume that {ug"', y^+i, Pg, w^j^^i solves 
and (P3n)) - (P3T1) forO<g<p<^-l. Then it is clear that we can immediately use 
1(3:46)1 - dsn]), and ((330)l -(jm]) to determine y^^i^, w^+i, and pp+i from {^^^^ ^^i, /5q, w^j^^i. 

We then substitute these into p.45p to determine u^'^^-y which completes the induction step. With the 
{u]^ TVq^iT PqTW^qYq^i determined, similar arguments show that we can use these along with the initial 

data {ull^.ull^} ^° ^^'^ """+2' P«+2; "^^f+ii ■^'f+2}- ^ 

Having constructed a solution to S^je^o — 0; we will use the implicit function theorem to find solutions 
to = for e > 0. However, to apply the implicit function theorem, we must first establish that the 
partial derivative of with respect to rji is an isomorphism. 
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Lemma 3.4. Suppose de e s t' G ^^af' '^'^'^ ^'^^ "^^J' defined by S6/f(-) = S(0, 6*^,-). 

r/ien the derivative 

is a linear isomorphism. 
Proof. Fix 

a, = (tA,°°2,V?^2,^£",^o",'^i+i,'^m,^o^i'S^<2) e Z;, 

and let 

Then from (|3.27[) - (|3.37[) . it is not difficult to see that the equation 

DEg^{r]£) ■ St]£ = ae 

is equivalent to the following system: 

6Xo = v"o' 



^0 > 
^0 — "^0 J 

„--,0J _ „;.00 



i'p+2 



C) 


■3d«+iln(a) 




2 dt-J+i 



5wo = cJq, 



Sii"'^ -I 



t=0 



and 



5pp+i - Qp{5pp, Swl^, (5wp) = Wp+i, 

where Ap'^, 0^, and Op are linear maps that depend implicitly on rjg. This system has the same structure 
as the system p.4ip - p.5ip from Lemma r3.3[ and a slight variation of the arguments used in the proof 
of the Lemma can be used to establish the existence of a unique solution for the given ag. In particular, 
this shows that D'E.g^ {r]i) is an isomorphism. □ 



By Proposition 13.11 the map : (— eo,eo) x g ^ x 3^^ ^-^ Zf is well defined and analytic for 
eo > small enough. Lemmas 13.31 and 13.41 then allow us to apply an analytic version of the implicit 
function theorem (see [6], Theorem 15.3) to conclude the existence of (shrinking eq if necessary) an open 
neighborhood (— eo,eo) x Uq^ C M x Xj^^^ of (0,0^) and an analytic map -qg : (— eo,eo) x Uq^ ~^ yjise 
that satisfies (e, 6e, f]e{£, — for all (e, 9i) e (— eq, eo) ^ Ug . □ 
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4. Post-Newtonian expansions 



4.1. The limit equation. Before discussing the e \ limit equation for the system (|2.62p . we first 
consider local existence and uniqueness of solutions to the cosmological Poisson-Euler equations. 

Proposition 4.1. Suppose s > 3/2 + 3 + £, S = ^io/{2Cs), Po £ BsiH"), and Wq G H^R^). Then there 
exists a maximal time Tq and a unique solution 

Pit) e Xt„ 

to the Poisson-Euler p.l0p - p.l2p with initial data p(0) — po + p and = /{po + p). Moreover, this 
solution satisfies 

{i\kt)w'{t))L^ ^ {nkmm)L^^ = 

for allte [0,To). 

Proof. As the system Poisson-Euler (|1.10l) - (|1.12p is clearly a (nonlocal) symmetric hyperbolic system, 
the statements concerning existence and uniqueness follow immediately from standard theory. To prove 

the second statement, we observe that ^pw^ satisfies 

dtipw-^) = -di{w^pw') — - d-' f{p) + pd-'<^ H ^pw^ 



a{t) 

Taking the inner product of this equation with 1 yields 



(4.1) 



1 /^/ttO| oJ.S.\ , 3a'(t) 0/, 



a{t) 



By prT2)) . we have a{t) {np\d'' ^) = {A'i>\d-' ^) and hence 

ait){np\d-'i>)L2 = - / djdi^d^^d^x = / dj{dii>d^^) d^x = 0. 
Substituting this into (|4.ip gives 
(4.2) dt{l\'f^')L2 



3a!{f) ooj 
= - ^ {l\pw') 



2a{t) 



L2. 



By assumption, (1| ipvo) |(=o)l2 = {\\'w^)i^2 — which combined with the differential equation (j4.2p shows 



that {l\l{t)w{t))L2 = for aU t G [0,To) 



□ 



From [27], we know that the appropriate e \ limit equation for the system (|2.62p is 

(4.3) AldtW ^ AidiW + Fa + C'dioj, 

(4.4) C'diW = 0, 
where 

(a{t) 0\ 
5^'^ 
10" 
\ 1/ 

The relationship between the Poisson-Euler equations (|1.10p - ()1.12p and the limit equation (|4.3p - (|4.4p is 
given by the following Proposition. Here and in the following section, we require the following evolution 



(4.5) 



^0 











and APq q 
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spaces 

i+i 

Proposition 4.2. Suppose s > 3/2 + 3 + 1, S = no/[2Cs), po 6 BsiH"), w^o ^ H^iR^), and let{p,w',i>} 
and Tq be as in proposition 14.11 Then 

W = {0,0,0,0, p,iO,w')) eXTo,e+2,s and = (,5^5^at$, 0, 0, 0, 0, 0) S A'to,^+3,s+i 
solve the limit equation (|4.3|) - (|4.4p . 





Proof. The proof follows directly from substituting W and ui into (|4.3p - (|4.4p while using p.l0p - p.l2p . 

^M-UM, mB-^M, and (i^TDD-dmn. □ 

4.2. Proof of Theorem 11.11 We are now ready to prove Theorem ll.il and thus establish the existence 
of cosmological post-Newtonian expansions to arbitrary order. 



Proof of_Theorem^ Given y^^ £ S4, y^' e §3, Po e Br{H^), wl 6 H^M?), iij-^ e H'^+^-'^iSs), and 
uli, e i?^-^(§3), let 

d?l^,u'/ = yl-' + eX-'{e), dfl^^u'^/ = e{Sy,J + u'^/{e)), 9f j^^^^^o = + = 1 

Pe|t=o = A*o + Ao(e) + Po, and (pw )|t=o - 



Aio + Ao(e) + Po 



be the initial data from Theorem 13.21 By construction, this data solves the constraint equations (|3.ip - 
p.3p . depends analytically on e, and satisfies (|3.13p - (|3.15p . In particular, this implies by (|2.5ip . (|2.52p . 
and Lemma EH that 



(4.6) We(t) = iulit),Wrjt),unt),cj^l^it),w,{t))^ 
satisfies 

(4.7) WM^C^ii-^o,eo),nn, 
and 

(4.8) Pf\t^o^e\\^.-,<l forp=l,2,...,£+l. 

Together, Proposition 14. 21 and (|4.6p - (|4.7p allow us to apply to Theorem 3 of [27] and conclude (shrinking 
eo if necessary) that for any T < Tq there exists maps 

We e XT.e+2.s < e < eo, 
p 

W e Xtj+2-p,s-p p^l,2,...,e, 
p 

We e XT,i,s-e-i {p, e) e ^>i+i X (0, eo), 

such that 

(i) We{t, x^) solves equation ([^IB^ on the spacetime region {t ^ ,x^) e M = [0,T) xT^ , 
p 

(ii) W (1 < p < t) satisfies a linear (nonlocal) symmetric hyperbolic system that depends only on 

q 
{W \ q — 0,1, . . . ,p ~ 1} where is a defined in Proposition [421 
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P 

(iii) We (p > ^ + 1) and satisfy the estimates 

\\W,{t)\\n^, + e\\dtWS)\\H^-^ < 1, and + e\\dtW,{t)\\^.-e-2 < 1 

for all {t, e) G [0, T) x (0, eo), and 

(iv) We admits a convergent expansion (uniform for < e < eq) of the form 

' P ^ P 

We^W + 22'^^'^ + iJ 

P=l p=£+l 

where the expansion is convergent in C"([0; T), W^^^'^). 
Finally, similar arguments as in the proof of Proposition 6.1 in [20] can be used to show that {ul^ = 
e~^ul^ , Pe,wl} determine, via formulas (|1.2ip - (|1.23p . a solution to the Einstein-Euler equations in the 
harmonic gauge and moreover that dtul^ = e'-^u^ ^ and diul^ = Wj ^ + dj^f . This combined with the 
statements (i)-(iv) above completes the proof. □ 
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